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Generalized Quantum Tunneling Effect and 
Ultimate Equations for Switching Time and 
Cell to Cell Power Dissipation Approximation 

in QCA Devices 


( N 
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owadays quantum-dot cellular automata has been playing a major challenge in competition of developing molecular 
■ ■ 

I circuits which overcomes certain fatal complications of CMOS devices for instance impurity inconsistency, high device latency - 

etc. Quantum-dot cellular automata empires the quantum binary logic gate designing technology by the help of fastest I 
U quantum bits. The electronic orientation assigns the binary q-bits inside the q-cells. In this letter one novel methodology for u 
.tunneling barrier resistivity analysis and device latency calculation have proposed. The tunneling barrier resistivity has . 
[propounded with a generalized mathematical expression to reveal the instantaneous effect of barrier resistance to the electron I 
j during tunneling. Proposed technique for switching time computation explores the signal propagation delay by considering u 

I the directional flow of quantum signals through different cells under different clock zones. Furthermore, in this letter one ■ 
novel cell wise energy and power dissipation computation technique has proposed which offers a very accurate power drop \ 
.data for the QCA wires. The proposed mathematical expressions are flexible enough for computing the power drop for any m 
w It onto axial wire and for any off centered or heteroaxial wire. This contribution extends a previous work on switching time, i 
[provided in reference. v v . , ' ^ ' l ? F? 'm. J 


Relevant Keys—Quantum computing ~ tunneling equations ~ switching time ~ device latency ~ poptr 

ou&p a c$nN> dissipation. ^ Wiip Wiip 


i. Introduction 

Quantum-dot cellular automata introduces the quantum 
logic gates with a controllable electron tunneling facility for 
the purpose of reconfiguring the electronic positions to attain 
the lightning device speed. The binary logic states are defined 
by the electron configuration change. The right handed 
electronic configuration terms the high logic state 4 1 ’ and the 
left handed configuration states the lower logic state ‘0 \fig. la 
[l]-[5]. fig. lb shows the QCA three input majority logic 
system, fig. lc is a clocked binary wire and the fig. ldisl cells 
inverter layout. 

In this contribution the tunneling barrier resistivity has 
mathematically calculated and formulized into a mathematical 
form which is able to compute the instantaneous behavior of 


the electrons during tunneling [6]-[10]. Furthermore, the signal 
propagation delay computation methodology has proposed 
which explores the delay of the signals in propagation from 
one cell to another. In this letter a novel power drop 
computation methodology has proposed which is able to 
monitor the power dissipation of each cell of a binary wire 
weather it is homoaxial or heteroaxial. This paper is an 
extended version of a previous publication on QCA device cell 
wise switching time approximation [2]. 


_ 

Author’s signature 

This document is self-published. The copyright restricts the commercial uses of this article but allows the citation. The author 
declares that he has full right to share the document publicly, this document is fully original and it has neither published before 

nor submitted at anywhere. 


S 


20 M H>rcor>n > no 












© THESSR DIGITALLIBRARY 

Sponsored by Physics Tomorrow 


Self-Publishing 


Page | 3 

T 
H 
E 

s 
s 

R 
D 
I 

G 
I 
T 
A 
L 
L 
I 

B 
R 
A 
R 
Y 


Q_ 


B Y=AB+BC+CA 


(a) Positive Negative 
polarization effect 


(b) 3 input majoriy logic 

_<§)|g)_<§)J1(§>_<§)J 

(c) A clocked binary wire 



(d) 7 cells inverter 


Fig. 1 (a) Positive/Negative polarization effect, (b) 3 input Majority gate, 
(c) A clocked binary wire, (d) 7 cells inverter 

ii. Aquired Contents From Previous Paper 

Layout of 2 bit binary to gray code converter. This is a 2 
bit code converter has designed followed by the layered 
design concept. Each atomic layer has separated by 11.50 nm 
from another. In this layout three layers have used for 
avoiding the wire crossover complications [2], [11]-[21]. 



Fig. 2 Two bit binary to gray code converter [2] 

This is the equation which states the tunneling rate in terms 
of frequency. The tunneling rate equation is given by [2], 

„ l 




Hz 


< 1 > 


This is S-instantaneous tunneling density (SITD), which 
defines the tunneling density at every instant inside the 
tunneling isolation barrier, 


lu = 


dr 
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dEg 



E u ^ u 
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hi. Analysis Of Tunneling Phenomena And 
Barrier Resistivity Calculation 

A. Calculation of barrier resistivity through tunneling rate 
and effective tunneling distance 

Below fig. 3 is the model of the tunneling barrier where V 
stands for barrier potential (in height) and r is the barrier 
thickness. 



The eq. 1 provides the numerical identification of the 
tunneling rate [2], which is inversely proportional to the 
tunneling resistance R T . The value of R T is dependent on the 
thickness of tunnel barrier 4 r’ and the potential height of the 
barrier ‘V’, fig. 3. The resistivity of this barrier is assumed as 
Pr, therefore the R T = ^- 7 =—. By replacing this value of R T 
in eq. 1, it is formularized as, 


r= - 
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is assumed as a constant value then the tunneling rate 


is dependent only on the variation of the 
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. Therefore, this is mathematically expressed 


which defines the proportionality. 

The width of the tunnel barrier is relative in nature. After 
starting the tunneling the effective value of r starts to decrease 
likewise r, r-Ar, r- 2 Ar, r- 3 Ar and so on, fig. 3. Along with the 
decrement of r the effective resistivity gradually reduces. 
Therefore, there are two different parts of eq. 2, one is a 
constant part and the second one is a variable part. 


T= ■ 


Sr***? t 


(1 - .i#) 
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r 0 is named as S-factor which is equals to 

This equation is re-expressed as where d is the effective 
tunneling distance [2], 

^(iT 

d L =r 0 c 1 
d=r 0 -1 c 

After replacing by the r D the equation becomes, d z 
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B. Solution for tunneling barrier resistivity pj 

Equation 3 is the simplified form of tunneling rate T. Let 
Page | 4 assume the barrier thickness is x nm. Hence, the effective 
tunneling path for the electron is x nm. The electron starts its 
journey from 0 nm and it finishes the journey at x nm, figJ. 
After integrating the eq. 3 in both sides with dx the resulting 
equation provides the speed of tunneling. The calculated 
tunneling speed of the electron is mathematically expressed 
as, 
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The quantum tunneling occurs at the speed of light. Therefore, 
in above equation, v=c= 3x 10 s m/sec. After replacing the 
value of v from above equation now it is rearranged as, 




Practically the value of r is equals to x as fig. 3, therefore, this 
equation is further extended as, 


Pr= — (l - e* 7 ) 


In purpose of point to point investigation of tunneling barrier 
resistivity r can be written as (r-nAr), where the effective value 
of the r is undertaken only. Ar is the small increment in 
effective tunneling distance. Assuming that Ar is a uniform 
increment. Moreover, n is the integer number which defines 
the current state (current position of electron inside tunnel 
barrier) of tunneling. It gives the values like 0, 1, 2, 3 ,, \ If the 
electron is its initial state (before tunneling) then the n=0, if 
next to initial state inside the tunneling barrier then n=l and so 
on. Therefore, 


Pr = 


[T-nAF\ 2 ( . 
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is found as, — 

rd 


From paper [2] the value of 

Or-P^l 

Equation 4 is the generalized form of tunneling barrier 
resistivity. This eq. 4 is named as S-tunneling resistivity. 
CASE 1: When n=0; 

Q m 


|(r) 2 rd| 

= 1—1 

1 sac 2 1 

\sac 2 1 


Pt MAX~ 

At this situation the electron faces the maximum resistivity 
from the tunnel barrier. 

CASE 2: When the electron is absolutely at the center of the 
tunneling barrier then the condition is written as, nAr=-^. 


_| W M |_| d r | 

MID I 2 I l-lEtf-C’ 2 ! 


Pt h 
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Therefore, at the middle of the barrier the resistivity is four 
times lesser than the initial resistivity. 

CASE 3: After the electron reaches to another dangling bond 
it is mathematically expressed as, nAr=r. 

Pr rest — ^ £Tm 

Tunneling barrier resistance R T is defined as, 

Here A is the effective area which is defined by rx V, as 

fig. 3. 

Therefore, R T =—; r defines the effective tunneling distance 

which is equals to d. This d is the width of the tunnel barrier. 
The tunneling occurs at speed of light. So, r is the product of c 
and t. Therefore this tunneling resistance is redefined as, 

c is the speed of light and t is the time taken by the 
electron to penetrate the barrier completely. By putting the 
value of^r Mjljr in^j equation, the maximum value of R T is 
expressed as, 

R , 




ct 

Hr 


Similarly, the intermediate tunneling barrier resistance is 
expressed as, 

□ = fL 

7 MID 

Now, the value tunneling resistance is put in eq. 1, 




r A(i 


- 1/2 


In a situation if in case, =0 then only the inverse of 
switching time will be the tunneling rate as following 
expression, 

r = - 

t 

From the calculations of section IV, the tunneling time is 
equals to the t in above equation. Therefore, the above 
tunneling rate equation is reformulated as, 

T = 5 X IQ 16 [l - ) Hz 

For QCA devices 50 X 10 1- Hz is a constant value which is 
named S-frequency. That is the QCA device operating 
frequency, from section IV. 

After a hard calculation the SITD has redefined as, 

E L.l 

, LC9^1D 3S .f 2KT 

lu=- 

T-T-sIkt 
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Fig. 5 Effect in DB under P+ (positive) polarization 


Alignment of electron 


0 1 2 3 4 5 6 

Tunneling barrier thickness or effective tunneling distance (nm) 

Fig. 4 Tunneling barrier resistivity analysis for several semiconductors by the help of above case studied equations. The pj . 

illustrated in above graph (according to case 1, case 2 and case 3 equations of barrier resistivity) 

This theory imparts, the Silicon has a barrier resistivity Pr MAX 
= 7.373624 Xl0 = Q cm. In standard, the resistivity of an 
intrinsic Si at 300K is 3.5 X 10" Q cm. This calculation 
proves that the resistivity becomes approximately double 
when the bulk material changes to respective the 
nanomaterial. Typically the resistivity increases when the 
temperature decreases as a result of decrease of the random 
lattice phonon. QCA devices can operate at maximum 17K 
temperature that avoids the electron ejection form the DB. 

Thereby, the resistivity increment is theoretically proven as a 
highly logical phenomena. 

iv. Scaling Of QCA Cells In Nanoscale 

In precise atomic view the Q-cells contain four q-dots 
which are typically formed by the dangling bond creation on 
the atoms are present in the lattice surface. In practice the 
removal of four hydrogen (H) atoms from four silicon (Si) 
atoms in (100) Hydrogen doped Si quantum dot results this 
dangling bond (DB) which are suitable for localizing the 
electrons [3], fig. 5, 6. Usually q-dots can have maximum 8 
nm to 9 nm diameter which are capable to localize 10 to 15 
atoms in it. But each Si atom can perform as a distinct 
quantum dot, this is the most versatile property of Si qdots, 
just because of that Si is highly preferable for QCA device 
fabrication purpose. But in QCA technology the q-dots have 5 
nm diameter [3], [20], that is the assembly of a very few Si 
atoms approx. 2 to 4. 


Pr MJJt and & are briefly 



H atom ^. 

Si atom 


Fig. 6 Effect in DB under P- (negative) polarization 


Under the effect of positive or negative polarizations the 
alignment of the electrons are shown in atomic view, fig. 5, 6. 
The blue circles represent the silicon (Si) atoms and the greens 
are hydrogen (H) atoms inside a quantum cell. Each corner of 
this cell has one quantum dot. Therefore, every cell contains 
four quantum dots. In the second row of this lattice two 
hydrogen atoms have removed from the 4 th and 5 th silicon 
atoms and similarly in the third row of this the lattice from 3 rd 
and 4 th silicon atoms two hydrogen atoms have removed, fig. 
5, 6. Due to the removal of 4 H atoms the four dangling bonds 
have created. The DBs are the unsaturated valance which are 
capable to localize the electrons. Using Scanning Tunneling 
Microscopy (STM) the precise electron implantation is 
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possible inside the DB. Those implanted electrons become 
free but immobile in natural condition (depolarized state). 
When the external polarization (in terms of electric field) is 
applied then the electrons gain sufficient energy for being 
tunnelable through tunneling barrier. 

Dimensional scaling of quantum cells 

One q-cell has 18x18 nm 2 area, fig. 7. Therefore, the 
length and the breadth both are equals to 18 nm. Out of 18 nm, 
(5+5) 10 nm space has covered by two upper q-dots. 
Remaining 8 nm space has a distribution like (1+6+1) nm. 
Each q-dot has the 1 nm separation from the end point (grain 
boundary), fig , 7. Therefore, the effective space between two 
q-dots are equals to 6 nm. 



Fig. 7 Scaling of the QCA cells in nanoscale 


The separation between two q-cells are 2 nm. Hence, the 
effective intermediate spacing between two q-dots of two 
different q-cells are given by 4 nm, fig. 7 [21]. 

The quantum tunneling occurs at the speed of light i.e. 300, 
000, 000 m/sec. Therefore, the tunneling time is calculated as 

t 1 = x ^ x I = 2x 10attosecond. 

The quantum signal propagation has the lightning speed 
during travelling from one cell to another cell. Hence, the 
inter-cell switching time is expressed as, £[ = 

£4 x 10 ^ sec ^ i 33 x 10 -z attosecond. 

The net switching time for two cells is formularized as, 
t s =2 t L +£i 

= 5.33 x 10 attosecond. 

F= — = 50 Peta Hz which is the operating frequency of the 

four q-dot QCA devices. One QCA cell is capable to be 
operated at 50 PHz frequency that is significantly higher than 
CMOS, NMOS and other electron island devices. 

v. Proposed Latency Diagram for Device 
Latency Computation 

In this section the latency of the proposed 2 bit binary to 
gray code converter has analyzed by the help of proposed 
latency diagram approach. This latency diagram is calculating 


the device latency theoretically with the proposed algorithm. 
In this process the whole circuit is subdivided into a few 
smaller segments where the signal flow is presented explicitly, 
fig. 8, 9. The general equation t s = (C — l)r _1 [2] for 
switching time computation have used to approximate the net 
switching time of this 2 bit binary to gray code converter 
device. Below diagrams show the calculation procedure of the 
device latency. 

Proposed algorithm for signal propagation delay (device 
latency) computation: 

Begin 

1. Select the inputs and outputs of the respective circuit. 

2. Draw the proper signal flow directions on the layout 
from inputs towards outputs. 

3. Group the cells based on the clock and their 
positions. 

4. Assign the weights of each group based on number of 
cells present in a group. 

5. Use this equation t s =(C — l)r _1 and calculate the 
switching time for each group of cells using the 
obtained weights i.e. C. 

6. Perform the simple addition of all the obtained 
results. 

End. 

After designing the cell to cell interaction phenomena, the 
following diagram is achieved which explains the signal flow 
direction within the q-cells, fig. 8. 



Fig. 8 Signal flow directions in 2 bit binary to gray code converter 
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Fig. 9 Device latency diagram for proposed 2 bit binary to gray code 
converter 

In above diagram fig. 9 , the code converter layout has 
subdivided into several group of cells depending upon the 
signal flow strategy which have undertaken to approximate the 
signal propagation latency during the device operation. This 


Self-Publishing 

diagram has some weighted values which is determined by the 
number of cells are present in a group as indicated above. The 
entire mathematical process have illustrated below. 

The general equation for switching time is t B =l(c- l^r- 1 1 
[2], Therefore, only with the known values of C and T this 
equation is solvable. From section IV, the value of T is 50 Peta 
Hz. The proposed latency matrix is given as 
L=[2 7 2 2 6 3 3 3] L;: .. ;g . which is generated from fig. 9. This 
matrix elements depend on the weights of this device latency 
diagram. 

Calculation 

t s 2 = (2-l)x 50 _1 = 0.02 attosecond. 
t s L0 = (10-l)x 50 _1 = 0.18 attosecond. 
t s L1 = (ll-l)x5Q -1 = 0.20 attosecond. 
t s 5 = (5-l)x 50 _1 = 0.08 attosecond. 

Hence the summed value of all of these provide the device 
latency or signal propagation delay. The sum of 
(0.02+0.18+0.20+0.08) is equals to 0.48 attosecond. 
Therefore, the value of t s is defined as 0.48 attosecond. It is 
established that this device has 0.48 attosecond signal 
propagation delay. 



(a) [9] 


(b) [ 11 ] 



(c) [12] 


Fig. 10 Latency diagrams for erstwhile designs by proposed methodology, (a) [9], (b) [11], (c) [12], (d) [13] 
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In above diagrams the device latency of the erstwhile designs 
have computed by proposed methodology. Fig. 10 a,b are the 
2 bit binary to gray code converter designs and fig . 10 c,d are 
the XOR gate designs. 
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and axial angular displacement between two cells, fig. 12. The 
kink energy of the cells increases if the liner cell displacement 
and translation angle 6 increases. This relationship is 
mathematically expressed as, 

E diss = hf 00 5in(n0) 
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Fig. 11 Device latency comparisons for previous designs (values from fig. 10 
latency diagrams) 

vi. Proposed Energy Dissipation 
Computation 

In QCA devices the energy dissipation is mainly dependent 
upon two major things, one is displacement between two cells 
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In this case the linear displacements among the cells have 
been illustrated by sketch. In above figure, r is the distance 
between two cells and gradually r increases with the 
increment of number of cells. Then it is denoted as, r'=(C-l)r; 
where C is the number of cells. From fig. 7, the computed 
value of r is 20 nm. 

The saved energy E s is expressed as, 

-Ea 


Fig. 12 Linear scaling of a QCA binary wire 

^=y(V- !} + 


sin(4?)\ 


( 6 ) 


E, = E T - L ((C - 1) + E; 


( 5 ) 


E T is the total energy that is provided for polarizing the driver 
cell. The saved energy is the amount energy that will be 
utilized to latch the next cells. 

Dissipate energy is given as, 


sin(as)/ 

In above equation / is the dimension of a single qcell and 
?" r =(C-l)r, is the separation between qcells. 

Case 1: For a homoaxial binary wire, the value of 0 
approaches zero that modifies the above eq. 5 as, 

E s =E 7 - f (C-l)£- a 

In above fig-12, if the cell 5 is selected which saved energy to 
be calculated then E s is written as, 

E s =E t -^CS-DF* 

E q = E t — 4.44E-. 
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Correspondingly, this approach is applicable for any cell that 
is present in any homoaxial binary wire. If any cell is 
displaced upward or downward then accordingly the value of 
r and 0 changes. Through this mathematical expression the 
power drop is also become computable. 

The dissipated power 

Pdiss is given as the ratio of the energy drop and switching 
time. This ration directly explicates the cell wise power drops. 



For a homoaxial wire, as fig. 12, the power dissipation is 
expresses by, 

i_!=" 

n _ L ^ _ J r 

— [-_i — 

For 18 X 10-mrr dimensional qcell the power dissipation 
equation is generalized as, 

^iss=2.78x 10 :3 r'E k W 

Assuming a standard value of the kink energy for a cell, 

£jt=3.14x 10 -E P J 

^*= 8722 . 22 ^ 

The above equation is independent to cell count, but it has a 
factor of distance (r) which is a function of number of linearly 
used cells. 

Case 1: For the second cell the value of r is expressed as r, 
as fig. 12. 

Then, 

F^,=8722.22r 

=8722.22x20x 10 " 5 
=1.744 mW 

Case 2: For the 6 th cell value of r is 5r, as fig. 12. 

^tsi=5rx 8722.22 
=8.722 mW 

Similarly, by this approach for other cells it is possible to 
calculate the dissipated power. 

Therefore, in case of a homoaxial wire, the power dissipation 
becomes independent of the number of cells. This only varies 
with the change in the value of r. If the cell count increases 
then the logical increment of the value of r comes about that is 
represented by r , in table I. 

In below figure one off centered QCA cell has shown. 



Fig. 13 An off centered/ heteroaxial cell 



TABLE I. Energy and Power Dissipation data for 6 th cell of 

THE QCA WIRE 


Sr. 

No 

Translation 
angle (S ) 

*"'=5r 

(nm) 

(J) 

J’d-. 

(mW) 

1 

0 

20 X 5=100 

- 

- 

2 

0.1 

100.5 

0.9642X 10” 3¥ 

9.642 

3 

0.2 

102.05 

0.9791 X 10" 3¥ 

9.791 

4 

0.3 

104.67 

1.004X 10" :¥ 

10.04 

5 

0.4 

108.57 

1.04 IX 10" 3¥ 

10.41 

6 

0.5 

113.94 

1.093 X 10" 3¥ 

10.93 

7 

0.6 

121.16 

1.162X 10“" 

11.62 

8 

0.7 

130.75 

1.254X UT“ 

12.54 

9 

0.8 

143.53 

1.377X ID” 3 * 

13.77 

10 

0.9 

160.87 

1.543X ID -35 

15.43 


Above table interprets the energy and power dissipations of 
the 6 th cell in the QCA wire due to the variation in translation 
angle (0" ), fig. 12. Through this approach the energy 
dissipation of any particular QCA cell can be achieved and the 
entire behavior of the wire can be justified mathematically. It 
is recommended that by keeping the notice on the dissipated 
energy to the last cell of the wire through this method, the 
driving energy should be provided to the driver cell. This 
could provide an idea of the effective power utilization and the 



Fig. 14 Effect of translation angle on power dissipation and energy dissipation 
of the 6 th cell of the QCA binary wire, fig. 12. 


TABLE II. Percentage Change in Power Drop for the variation 

of Translation Angle 


Sr. No. 


Translation angle 

(O 


Percentage change (%) 
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Above table II, enlightens the change in power dissipation in 
percentage for the 6 th cell if the translation angle variation 
occurs, remarkably. This analysis justifies the gradual 
increment of the dissipated energy compared to the homoaxial 
wire, due to the increment in the translation angle. This 
section concludes and proves the power dissipation 
characteristics of a QCA cell in a wire and justifies the 
deviation in power drop for an off centered cellular wire. 

vii. Conclusion 

In present quantum nanotechnology there are many 
instances alike Carbon Nanotube FET (CNTFET), SpinFET, 
Single electron transistors (SET), Electron Spin Devices etc. 
which are highly adaptive in performance. Quantum-dot 
Cellular Automata is also a molecular logic synthesis 
technique which is fastest in performance than others. The 
quantum tunneling effect gives it the lightning speed of 
operation. In this research letter two major things of QCA 
devices have analyzed. One is the tunneling barrier resistivity 
and another one is the switching time. It is mathematically 
proved that the variation of the barrier resistivity causes the 
variation in electron tunneling probability which affects the 
tunneling rate notably. The proposed S-tunneling resistivity 
equation is enough to explain the instantaneous behavior of 
the tunnel barrier during tunneling. The novel switching time 
commutation technique provides a major sense in QCA 
devices by calculating the switching time accurately. 

In section III, the tunneling barrier resistivity equation has 
proposed which reveals the nature of electron tunneling during 
the penetration of the electron through the barrier. In the next 
context the switching time computation procedure has 
proposed through a mathematical model which explicates the 
signal propagation delay along the Q-cells. This signal 
propagation time is mathematically indicated by t s which is 
approximated by several complex calculations in section IV 
and V. Proposed theory imparts the operating frequency of the 
device is 50 Peta Hz and the time required for an electron to 
tunnel the barrier is 2xl0 -2 attosecond. By the help of 
proposed latency computation methodology the device latency 
of 2 bit binary to gray code converter has solved and the 
latencies for previous designs explored from erstwhile articles 
have also illustrated in section V. This article contributes more 


by presenting a novel energy and power dissipation analysis 
methodology. In section VI, the dissipated power equation has 
proposed that able to calculate the change in power drop due 
to the angle variation in case of off centered cell. It is reported 
that if 0 . 1 ' angular translation occurs then the power drop 
rises to 10.54% and so on for the further angle increment. 
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